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A new geometrical perturbation scheme is developed in order to calculate the electromagnetic fields 
produced by charged sources in prescribed motion moving in a non-straight perfectly conduct- 
ing beam pipe. The pipe is regarded as a perturbed infinitely long hollow right-circular cylinder. 
The perturbation maintains the pipe's circular cross-section while deforming its axis into a planar 
qq , space-curve with, in general, non-constant curvature. Various charged source models are considered 

f"^ ■ including a charged bunch and an off-axis point particle. In the ultra-relativistic limit this permits 

a calculation of the longitudinal wake potential in terms of powers of the product of the pipe ra- 
£SJ , dius and the arbitrarily varying curvature of the axial space-curve. Analytic expressions to leading 

order are presented for beam pipes with piecewise defined constant curvature modelling pipes with 
straight segments linked by circular arcs of finite length. The language of differential forms is used 
throughout and to illustrate the power of this formalism a pedagogical introduction is developed 
, by deriving the theory ab-initio from Maxwell's equations expressed intrinsically as a differential 

' system on (Minkowski) spacetime. 

PACS numbers: 87.56.bd, 02.40.-k, 31.15.xp, 41.20.-q, 29.27.-a, 41.60.-m, 41.75.-i 
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Considerable activity is being devoted to the design of advanced machines that can produce pulsed sources of 
intense focussed electromagnetic radiation. Such sources offer unprecedented opportunities for probing the temporal 
q ■ and spatial microstructure of processes in Nature. Many designs rely on being able to control the motion of ultra- 
• ' relativistic electron beams by external fields in beam pipes with non-uniform spatial curvature. The production 
of femtosecond radiation pulses requires high peak electric currents and the maintenance of low emittance electron 
' beams for the self-amplified spontaneous emission of X-rays demands extreme design criteria in order to sustain beam 
O I stability in the presence of radiation back-reaction on accelerated sources. 

A direct analytic approach to this electrodynamic problem via the coupled system of Maxwell's field equations and 
the equations of motion for the particle beams encounters difficult problems due to non-linearities and retardation 
effects [H H, Q. A direct approach using statistical methods suffers from similar complications. Numerical approx- 
imations based on these equations exist but few are able to address the full complexities encountered in a realistic 
t^- , situation. In particular even when non-linearities are deemed negligible, the effects of (conducting) boundaries on the 
accelerating source via the radiation fields are often ignored in an attempt to build tractable models and the effects of 
pipe curvature are often restricted to those produced by motion in arcs of circles 0. Although insight can be gained 
OO ' from such modelling we feel that much remains to be understood in more general scenarios. 

In this paper we explore a new geometrical perturbation scheme that addresses some aspects of the general elec- 
■ trodynamical problem of charged sources in prescribed arbitrary motion moving in a non-straight beam pipe. In the 
analysis below the pipe will be regarded as a perturbed infinitely long hollow right-circular cylinder. The perturba- 
tion will maintain the pipe's circular cross-section while deforming its axis into a planar space-curve with, in general, 
non-constant curvature. It is assumed that the curvature k(z) depends on the arc-length \z\ of the planar space-curve 
and tends to zero as z — > ±oo. Furthermore if the cylinder has fixed radius a we require that |/t(z)a| -C 1 for the 
perturbation analysis to be effective. Under these conditions the initial-boundary value problem for the Maxwell fields 
in a perfectly conducting hollow perturbed cylinder given prescribed currents will be considered. The source currents 
will excite superpositions of electromagnetic modes of the empty perturbed cylinder as well as generating their own 
"self-fields" . The latter will include acceleration-fields induced by the motion of the sources in the curved regions of 
the beam pipe. A perturbation scheme will be established to calculate all fields in the pipe as a perturbation expansion 
in powers of an. Mode expansions based on Dirichclct and Neumann cigen-functions of the Laplacian for a circular 
disc domain can be used to reduce the general problem at each order to a two-dimensional linear telegraph-type 
equation with prescribed sources. The general causal solution of this equation is known and from it the complete set 
of electromagnetic fields can be constructed. Thus the perturbed beam pipe impedances can be found. Furthermore 
the fields arise naturally in a space-time description and thereby offer direct input into "leap-frog" coding schemes 
that couple the Maxwell sector to the equations of motion of the sources. 

In section 2 the notation used in the paper is established and illustrated by means of a pedagogic review of exterior 
methods used for solving Maxwell boundary value problems in sections 3 and 4. Section 5 deals with the introduction of 
Dirichclct and Neumann modes used to reduce the Maxwell system with perfectly conducting boundaries to a telegraph 
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type equation in two dimensions. Section 6 uses a geometric perturbation approach to explore the dependence of the 
radiated power from a smooth longitudinal convective current on local beam-pipe curvature. Section 7 deals with a 
moving point charge source and discusses in some detail radiant instantaneous power. This permits a perturbative 
calculation of the ultra-relativistic wake potentials in a planar pipe with arbitrary curvature and explicit analytic 
expressions can be found in the case where the curvature is piecewise constant modelling pipes with straight segments 
linked by circular arcs of finite length. Finally section 8 considers the radiation from smooth convected localised 
bunches with fixed total charge while the Appendices tabulate coupling coefficients and transfer kernels used in the 
main text. 



II. NOTATION 

The exterior calculus of differential forms offers a versatile and powerful means for analysing Maxwell's equations 
0, Q. The notation used below follows standard conventions for a manifold M with a metric tensor field. Thus TTM 
denotes the set of vector fields and TA P M the set of p-form fields on M. Metric duals with respect to any metric tensor 
g are written with a tilde so that X = g{X, -) G TA 1 M for X G TTM and a = g^ 1 (a, -) G TTM for a G rA 1 ^". 
The Hodge dual map associated with g is denoted by a star so that the canonical n-form measure ( "volume element" ) 
on an n-dimcnsional manifold M is the image of 1 under the Hodge map. In 4-dimcnsional spacetime we adopt the 
flat Minkowski metric. In a 3-dimensional space we adopt the Euclidean flat metric and regard time as an evolution 
parameter for forms in three dimensions. In a 2-dimensional space we adopt the Euclidean flat metric and regard time 
and a longitudinal coordinate as parameters for forms in two dimensions. One must then distinguish notationally 
between the different metrics introduced and their associated Hodge maps. However for any manifold M with Hodge 
map * one always has the standard relations 

$A** = *A*$, for <P,^eTA p M (1) 

«x*$ = *(*AX), for X G TTM, $ € TA P M (2) 

where ix denotes the interior (contraction) operator on forms. 

Maxwell's equations find their most cogent formulation as a theory of 2-forms on spacetime modelled on a space 
and time oriented 4-dimensional manifold with a metric tensor field g of Lorentzian signature (— , +, +, +). On such 
a spacetime M the set {e°, e 1 , e 2 , e 3 } will denote a local g-orthonormal coframe (a linearly independent collection of 
1-forms). The Hodge map associated with the Lorentzian metric g will be denoted by 7k-. Then 

« X $=- + $AI, for X G TTM, $ G TA P M (3) 

**$ = (-l) p+1 $, for $ G TA P M (4) 

For manifolds with a Euclidean signature and different dimensions these last two relations change as will be indicated 
for three and two dimensional spaces below. Finally note that for all n-dimensional manifolds of any signature one 
has the useful results: 

ix® A * = A ix®, for $ G I\A P M, 

* G rA 9 A/, p + q > n + 1 (5) 

d$A* = (-l) p+1 $ A + d($ A for $ G TA P M , 

* G rA ? M (6) 

As a notational convenience the expression below implies that "J" £ TA q N is a differential q-form on N where 

(<?) 

the manifold N follows from the context. 



III. ELECTROMAGNETIC FIELDS IN SPACETIME 



Maxwell's equations for an electromagnetic field in an arbitrary medium can be written 

dF = and d-kG = j, 



(7) 
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where F G TA 2 M is the Maxwell 2-form, G G TA 2 M is the excitation 2-form and j G TA 3 M is the 3-form electric 
1J| . To close this system, "electromagnetic constitutive relations" relating G and j to F are necessary. 



current source 



The electric 4-current j describes both (mobile) electric charge and effective (Ohmic) currents in a conducting 
medium. The electric field e G Th^M and magnetic induction field b G TA l M associated with F are defined with 
respect to an arbitrary unit future-pointing timelike 4-velocity vector field U G YTM by 



F 



and 



c b = iij-kF. 



Thus iue = and itjh 
Since g(U, U) = -1 



F = eAU -*(cobA£/). 



(8) 



(9) 



The field U may be used to describe an observer frame on spacetime and its integral curves model idealised observers. 
Likewise the displacement field d G TA X M and the magnetic field h G rA 1 ^/ associated with G are defined with 
respect to U by 



Thus 



d = ijjG , and h/co ~iu*G. 



G = d AU — * ((h/ Co ) A J7), 



(10) 



111) 



and i[/d = and zyh = 0. It may be assumed that a material medium has associated with it a future-pointing timelike 
unit vector field V which may be identified with the bulk 4-velocity field of the medium in spacetime. Integral curves 
of V define the avera ged world-lines of identifiable constituents of the medium. A comoving observer frame with 
4-velocity U will have [3] U — V. 



IV. TIME DEPENDENT MAXWELL SYSTEMS IN SPACE 

On any n-dimensional manifold a chart sets up a correspondence between points on some region (patch) on the 
manifold and a set on M. n . Thus in a 2-dimensional patch let £ = (^,^ 2 ) be a generic set of coordinates. Similarly 
let £ = ($ 1 ,^ 2 ,^ 3 ) denote coordinates on a patch of a 3-dimensional manifold and £ = denote coordinates on a 

patch of 4-dimcnsional spacetime. 

Let d denote exterior differentiation in any domain of a Euclidean space with coordinates £. Similarly let d denote 
exterior differentiation in a patch with coordinates £. A "moving" orthonormal (Cartan) coframe in flat spacetime 
with Minkowski metric g is a set of (independent [l6|]) 1-forms {e°, e 1 , e 2 , e 3 } with e° timelike. In general this will 
depend on the choice of coordinates £ in the sense that its exterior derivative will not be zero. In the following we 
adopt an incrtial frame with laboratory time £° = t and e° = cq dt with {e 1 , e 2 , e 3 } independent of t. Thus in general 
the coframe "moves" as a function of £. If [3 is any form on spacetime it will be convenient to adopt the abbreviation 
$ for Ca_f3, where Cx denotes the Lie derivative @, B, with respect to X. Thus e k = for k = 1,2,3. Within this 

dt 

framework introduce the tensor fields: 

g = e 1 <g> e 1 + e 2 ® e 2 , g = g + e 3 <E> e 3 , g = -e° ® e° + g, 

where e° = Co dt and <? is the metric tensor field on Minkowski spacetime. At each instant (t =constant), g is the 
induced metric tensor on Euclidean space and g is the induced metric tensor on the 2-dimensional submanifolds 

(leaves) where £3 = constant. Denote the Hodge map associated with g by # with 
and that associated with g by ff with 

#l = #lAe 3 . 

Then 

*1 = #1A e° := e 1 A e 2 A e 3 A e°. 
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To accommodate the effects of signature it is convenient to introduce the involution operator r\ on p-forms $ by 
?7$ = Then 

**=-V, ## = 1, ## = ??■ (12) 

By linearity the action of the Hodge map on an arbitrary form in Euclidean 3-space readily follows by expanding 
it in an orthonormal basis and using the relations 



#e 1 


= e 2 


Ae 3 




= e 3 


Ae 1 


#e 3 


= e 1 


Ae 2 



on the basis forms. Furthermore in a 2-dimensional Euclidean space 

If j3 (£) is a p-form on spacetime but generated by forms in the exterior algebra generated by {e 1 (^), e 2 (f ), e 3 (£)} 
(p) - - - 

then at any event with coordinates £ one has 

(p) V 

where, for each multi-index I, the set of exterior p-forms {e 1 (£)} denotes a basis for p-forms generated from the set 

{e 1 ^), e 2 (£), e 3 (£)}. One refers to the functions /?/ as the components of (3 in the e 1 basis. With this notation 
~ (p) 

(p) / ^ 

Define the 2 + 1 split of /3 (£) into the pair { /3 (£), /3 (£)} by the unique decomposition with respect to tf £ 3 : 
(p) (p-i) (p) 

0(0= /3 (e)Ade 3 + /3(o, (13) 
(p) (p-i) (p) 

where j3 (£) and /3 (£) are p — 1 and p- forms respectively, generated from the 1-forms in {d^ 1 , d£ 2 } satisfying 

(p-i) (p) 

ia /3 (£) = and i a /3 (£) = 0. Thus /3 and /3 are forms that do not contain d£ 3 . 
^(p-i) ^(p) (p-i) (p) 



It follows that for q = 0, 1, 2: 



#(/3 Ae 3 ) = #M), (14) 

(9) (?) 



#(^)=#(/3)Ae 3 . (15) 

(?) (9) 



For any 0-form (3 

(o) 



d P =d p +(£ a (3)dt z , 

(0) (0) ^(0) 
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where 

(o) o^ 1 (0) d£ 2 (0) 

From this it follows that, for q = 0, 1, 2: 

d /3 = d /? + d£ 3 A (C a 0), 

(?) («) W (<l) 

where d acts[l7} on exterior forms generated by {d^ 1 , d£ 2 }. Note that for all 2-forms one has d = 0. Let the 
3 + 1 split of the 4-current 3-form be 



(2) (2) 



with is J =0. Then, from Q 

8t (2) 



yields 



J (0 = -/,(0Adt+ p(0#l, (16) 

(3) (2) (0) 



dj = 0, (17) 



dJ(£)+ p(0#l = 0. 

(2) (0) 

It is convenient to introduce the (Hodge) dual forms: 

E:=#e, D:=#d, B:=#b, H:=#h, 

(2) (1) (2) (1) (2) (1) (2) (1) 

so that the 3 + 1 split of the spacetime covariant Maxwell equations (0) with respect to d t becomes 

d e = — B, (18) 

(1) (2) 

dB=0, (19) 

(2) 

d h = J + D, (20) 

(1) (2) (2) 

dD = p#l. (21) 

(2) (0) 

Allp-forms (p > 1) in these equations are independent of e° but may depend on t. Furthermore they are independent 
of the choice of (stationary) spatial co-frame constructed from {d^ 1 , d£ 2 , d£ 3 }, in any chart with local coordinates 

In the following it is assumed that b = u h and d = e e (with constant s, u ) where e — e r en , u = u r Mn ■ Thus 

(i) (i) (i) (i) 

in terms of e , h , E , H : 



(1) (1) (2) (2) 



d e = -uH, (22) 

(1) (2) 



dH = 0, (23) 

(2) 



d h = eE + J, (24) 

(1) (2) (2) 



ed E = p#l. 

(2) (0) 



(25) 



Electromagnetic Fields in a Curved Beam Pipe 



6 



V. THE MAXWELL SYSTEM WITH SOURCES IN A CURVED BEAM PIPE 

In terms of time dependent 1-form electromagnetic fields with general sources ( p , J ) the Maxwell's equations in 

(0) (2) 

a medium with scalar permeability p can be written 

d e + /!# h = 0, (26) 
(i) (i) 

d h pY 2 # e - J = 0, (27) 

(1) (1) (2) 

d#h=0, (28) 
/iF 2 d#e - p#l = 0, (29) 

(1) (0) 

where the admittance Y = 1/ (/zc) with c being the speed of light in the medium. 

These equations involve time dependent forms and are independent of particular local spatial coordinates. They 
depend explicitly on the Euclidean metric and for conserved sources define a well posed initial-boundary value problem. 
We now choose a coordinate system adapted to the interior Li of a beam pipe with a circular disc cross-section of 
fixed radius a at every point and an axis given by a planar space-curve with, in general, non-constant curvature k and 
\kcl\ -C 1. At each point on this curve one may erect a triad of orthogonal vectors in space, one member of which is 
tangent to the curve. The remaining vectors define a transverse plane. All points in the interior hi of the beam pipe 
lie on some transverse plane associated with such a triad with origin at some point on the axial space-curve. Let the 
region U CM 3 inside the beam pipe be described in terms of coordinates £ 2 , £ 3 ) := (r, 8, z) adapted to the central 
space-curve with curvature n(z) such that 

< r < a, < 8 < 2tt, -co < z < oo. 

A convenient field of orthonormal coframes Q on U is given in these coordinates by 

{e 1 = dr, e 2 = rd8, e 3 = (1 - eKo(z)x{) d z}, (30) 

with x\ = rcos8. Thus the Euclidean metric tensor g on hi is given by 

g = e 1 ® e 1 + e 2 ® e 2 + e 3 g) e 3 . 

In these coordinates the pipe boundary is the surface r = a, the coordinate z measures arc-length along the space-curve 
and on the space-curve r = 0. 

It proves convenient in the following to write k(z) = eno(z) and use e as a book-keeping device in order to keep 
track of different orders of k. In terms of adapted coordinates 

g = g + (1 — eno(z)xi) 2 d z <g> d z, 

where for each cross-section at constant z one has the induced metric tensor g on the 2-dimcnsional disc (0 < r < 
o,0 < 6 < 2tt): 



9 



dr®dr + r 2 d9 ®dQ. 



The associated contravariant tensors are 

-l 



d d Id d 
dr dr r 2 88 88 ' 



a- 1 = 9- 1 



1 8 8 



(1 — eKo(z)xi) 2 8z 8z 

The source forms will be expressed in terms of the scalar functions J r , </#, Jq and p of (e, t, z, r, 8). We choose to write 



J as 

(2) 



J (e,t,z,r, 8) 

(2) 



Jg(e, t,z,r,6)dr + rJ r (e, t,z,r,8)d8j Adz(l- eno(z)r cos 8) 
+Jo(e,t,z,r,8)rdr Ad8, (31) 
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from which it immediately follows that 

# J = —Jgrd9 + J r dr + J (l — eKo(z)xi)d z, 

(2) 

with orthogonal components {—Je, Jr, Jo) of # J . The associated electric current vector field is 4t J . 

(2) (2) 

The objective is to solve (l26l) . (1271) . (l28l) . (l29l) for the fields e and h on W in terms of prescribed sources and 

; (i) (i) 

initial data as a perturbative expansion in the axial curvature of the beam pipe. The strategy will be to project 
the field system into suitable modes that ensure that perfectly conducting boundary conditions are satisfied at the 
pipe boundary. In the adapted coordinate system this is achieved with the aid of complex Dirichelet and Neumann 
eigen-modes of the two-dimensional Laplacian associated with each transverse plane in the beam pipe. 



A. Dirichelet Modes 



Let T> be the smooth 2-dimensional submanifold (z=constant) with boundary dT>, embedded in Eucldean M 3 . The 
tensor g on D is that induced from the Euclidean metric g in M 3 . A complex Dirichelet mode set {$>n} is a collection 
of complex eigen 0-forms of the Laplacian operator d jfd on T> (associated with the metric g and Hodge operator 
#) that vanishes on &D . This boundary condition and the nature of the domain determine the associated (positive 
non-zero real) eigenvalues flfj. The label N here consists of an ordered pair of real numbers. Thus 

$ N :V^R, r,6^$ N (r,6), (32) 

satisfies 

d #d $jv + 0n®n#1 = 0, (33) 
with $7v|ax> = 0. It is straightforward to show from these properties that if 0fj ^ ^ 0, (N ^ M) then 

177*at#i = 0, 

where the bar denotes complex conjugation. If one normalises these modes so that 

*^*n#1=J*nSnm, (34) 



v 



V 



then it is also easy to show that 



An explicit form for <&jv is for neZ 



/ d $ N A #d $ M = 0nM% 5 Ni 
Jt> 



N (r,6) = J n (x q{n) -^e m \ (35) 



where J n {x) is the n-th Bessel function and the numbers {x q ^} are defined by Jn(%q(n)) = an d N := {n,q(n)}. 
The eigenvalues are given by {(3n = x q ^/a\. It follows from the integral flpj ] 



2 

drrj m (^ X q( m )~^) \^ 1 '^"^q) ~2^rn+l{ X q(m))^q(m),q'(m)^ 



that J\f% = 7ra 2 J^ +1 (x g( „)). 
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B. Neumann Modes 

In a similar manner one defines a Neumann mode set n} as a collection of eigen 0-forms of the Laplacian operator 
on D such that #e£ <3/jv vanishes on dT>. This alternative boundary condition and the nature of the domain determine 
the associated (positive non-zero real) eigenvalues a 2 N where again the label N consists of an ordered pair of real 
numbers. Thus 

satisfies 

d #d *jv + c4#iv#l = 0, (36) 
with #d ^ n\ov = It is straightforward to show from these properties that if ct 2 N ^ a 2 M ^ 0, (N ^ M) then 

^7* w #l = 0. (37) 

If one normalises these modes so that 



*m*jv#1 = M 2 N S NM , 



p 



then it is also easy to show that 



An explicit form for is for m £ Z 



d *at A #d * M = a 2 N M 2 N 5 NM - 



9 M (r,e) = J m .(x; (m) ^) e »» e , (38) 
where the numbers {%' p r m -\} are defined by J' m { x ' p ( m )) = ^ and M '■= {m,p(m)}. The eigenvalues are given by 

C. Mode Decompositions 

Since IA is simply connected one can represent the electromagnetic 1-forms e = e (e, t, z. r, 0) and h = 

(i) (i) V ' (i) 

h(e,t,z,r,6) as [n| 



e(e,t,z,r,6) = ^ V^(e, i, «)d$jv + X] V$(e,t,z)#(dz A 

2V AT 

+ ^ 7 ^(e,M) < Mr,0) t /z, (39) 



(i) 



TV 



N M 



+ ^2^(e,t,z)^ M {r,e)dz. (40) 

M 

Here, for any scalars /m with M = (m,p(m)) and Jn with TV = (n, q(n)), the summations above are abbreviations 
for: 

X f M " ' = X! X /m,p(m) ■ ■ • : X ^ " ' = X X -^>«(™) 
M meZp(m)6N N neZg(n)£N 

For future convenience the further abbreviation: 



X f N " ' = X X ^.?(«) 

JV n=m±lq(n)GN 



Electromagnetic Fields in a Curved Beam Pipe 



9 



will be used. The expansions above in terms of $at and ^ m and their derivatives ensure that the electromagnetic 
fields satisfy perfectly conducting boundary conditions at the surface r = a. 
Since $| r = a = we note that 



and for m, n G Z 



Furthermore with n G Z, 



d$ M Ad$ w = / $jvd*M = 0, (41) 



2tt 



f — 

Jo 2lT 



so 



Hence 



and 



J- n (x) = (-l) n J n (x), 
p(m) = p(— m), and g(n) = g(— n). 

^n,g(n) " n.q( — n) 5 $n,q{ri) ' $—n,q{ — ri)i (^^) 



d*n l9 (n) = $_„,,(_„), d* m , P ( m ) = (-l) TO d*_ OTiP( _ m) . (43) 

From Eq. (|43p and the reality conditions, 



e = e , and h = h , 

(i) (i) (i) (i) 



one has 



V E = (— l) n V E V H = (—l) m V H 

n,q(n) \ I — n,q(~ n) ' m,p(m) \ J — m,q( — m)> 

'n,<j(") V ^ r_„ jg (_ n )) 



(44) 



T E =( — '\\ n T E T H =(—'\Y n T H 

± n,q{n) \ x ) - t _„ )9 (_„); x m,p(m) V ^ - 1 — m,q{— m)' 

= f — lVV^ 

Similarly from the relations 

zJ^(z) = nJ n (z) - z.J n+1 {z), 

it follows that 

Mln.X-n) =K,q(nV and ^-m^-m) = ^^p(m)- ( 45 ) 

These relations enable one to pass simply from complex to real representations of the mode summations for the fields 
above. 

D. Perturbation Expansions 

Since for small |ko| the beam pipe approximates a straight cylinder we adopt the perturbative field-mode expansions 

V$(e, t, z) = V* i0) (t, z) + eV* (1 \t, z) + 0(e 2 ). (46) 
I§(e, t, z) = I^°\t, z) + el^\ t , z) + 0(e 2 ). (47) 
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T§ (e, t, z) = 7 ^ (0) (t, z) + 6 7 £ (1) (t, z) + 0(e 2 ), (48) 
with analogous expansions for the magnetic modes Vj^ , ImiIm anc ^ express the sources as a power series in e: 

J e (e, t, z, r, 9) = J^ 0) (t, z, r, 0) + eJ e (1) (t, z, r, 0) + 0(e 2 ), 
J r (e, t, z, r, 0) = (t, z, r, 0) + ejW (t, z, r, 0) + C(e 2 ), 
J (e, t, z, r, 0) = J (0) (i, z, r, 0) + eJ^t, z, r, 0) + 0(e 2 ), 
p(e, t, z, r, 0) = p(°) (t, z, r, 0) + ep«(i, z, r, 0) + 0(e 2 ). 

These expansions are then inserted into P?)) . ([%!?)) and 0(e n ) systems extracted for n = 0, 1. For general 
sources it is somewhat tedious to project out the equations for the perturbative field coefficients above. This is 
achieved using the orthogonality relations between the different Dirichelet and Neumann mode sets and the explicit 
relations listed in Appendix A to integrate (|26|) . ([27]) , (|28|) . (|2Q[) over the domain T>. 

In order to express the resulting equations for the mode amplitudes in a unified way we define 



'M AT {J~, Q, x, y) := J drr k T m [x p[m) ^ Q n (y q ( n ) ^) 

and note 



Furthermore if T is a Bcsscl function J or its derivative, J': 



^m,n(F,F, x, x) — ( l) m+n a_ m ^_ m ^_ n ^_ n - ) (J r , -T 7 , x,x). (50) 



Thus 



drrJ m (xp( m ) — J J' n ( x 'q(n)~) — ^M,n(J' J'i 



The symbol N (!F,G,x,y) enables one to write more compactly certain overlap coefficients that arise in the pro- 
jections of (|27[) . (|28p. (|2U)) over P. These coefficients are given in Appendix A. 
Using these definitions the projection of (|26|) to lowest order 0(e°) is 

^ (0) ' + M^ (0) -7^ (0) =0, (51) 
^'-^=0, (52) 
»iT - <^ (0) - 0, (53) 

and to first order 0(e 1 ) is 

AT AT 



/ T/ /f(l)/ r-H"(l)\ 2 A/(2 , / ^Y V A/ ff ( )' i f # (0) \ 77.$ 
(^Af ~ M J M + K o(z) 2^ 0V + /"V ) Af , AT 

Ar 



(55) 



iV 



2 T/ ff (°h Z?* 

Af.JV 

N 

wft (0) (^-<,) = o. 



(56) 



-«o 

N 
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where a superscript prime attached to a function denotes its partial derivative with respect to 
functions in ([M]) , (|55j) . ([55)1 depend on z or z,t. 
The projection of ([27]) to lowest order O(e ) is 

+ I jf\#dW^) Adr+ [ r4°\#dW^) A dO = 0, 

(4 (0) ' + ^^M (0) )/3^ 
+ / jf ] dW^Adr+ [ rJ^dW^ Ad9 = 0, 



and to first order 0(e 1 ) is 



J T> 



mi 



-{iT'-»y 2 < {1) -< {1) WmM* 



M 



N 



N 

+ K ' (z)Y^In {0) Gm% + I (4 X) " 4 O) K a (z)rco S 0)(#dW^) A d 



+ / r{4 1] - 4 0) K (z)r cos 6)(#d^ M ) Ad6 = 0, 
Jv 



N 

+k'o(z)Y'iT )f mn+ [ (J^ ~ 4 0) K Q (z)r cos 9)dW^ A dr 

N Jv 

+ I r(4 1] - 4 a) K {z)rcos9)dW^ Ad8 = 0, 
Jv 

N 

+Ko(z)J2' I N Q) (Cf LN ~Dt LN )- [ jW$77#l = 0. 

N Jv 



The projection of ([2"8)1 to lowest order O(e ) is 
and to first order 0(e 1 ), 



H(0)> 2 T H(0) n 



M.N 



N 



Finally the projection of ([25)1 to lowest order O(e ) is 
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and to first order 0(e 1 ), 

^ 2 (7m (1) '-^m (1) )^- / (p (1) -p (0) K o(^)rcos^#l 

Jv 

AT V > 



(66) 



1. Decoupling to lowest order 

The equations above to lowest order describe the fields that can be excited by sources in a perfectly conducting 
straight beam pipe. They are readily decoupled by substituting ([55)1 and (|6"5j) into ([57]) yielding an inhomogeneous 
telegraph-type equation for "f^^ ■ 

-H(0) 2 H(Q)n 2 2 H(0) 

In ~ c In + c a Nl N 

4°\#Wn~) ^ dr + J v r4°H#W N ~) A dflj. (67) 



X 2 



'iV L 

Similarly the equation for 7^ follows by substituting |59|) and ([65]) into ([5Tj) 



Once these equations are solved the other field components can be calculated from 7^ and 7^ and their 



derivatives. Thus ([65 ]) . ([53 ]) . <[59 ]) and (JG3J) yield 



VS m = 47j (0) , (70) 

_ 1 f,.v2-E(0) , 1 /* t(0) 



^ W = -jjr (^7* + J7 rJ v JT**r#iJ , (71) 

^ (0) = ^7^ (0) '- (72) 

a 7V 



These solutions are consistent for sources satisfying ([T7| . 
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2. Decoupling to First Order 

The first order equations involve the solutions to the lowest order fields. Similar to the case above the equation for 
7^ (1) follows by substituting (55]) and (JMJ) into pi)) : 



..H(l) 2 H(l)ti 2 2 H(l) 



C 2 K (z) 



Ml 



M.N 



M\ 



k' (z) V'/" b(0) (7*'* 

K 0^J 2^ JV M,N 



N 



M 2 



M JT> 
2 



M 2 



M JT> 



(4 1] - J { g 0) K {z)r cos6»)(#d* M ) A dr 
r(jW - J, (o) Ko (^cos60(#cZ*m") A<i(9 



iV 



A4 M 5z 
kq(z) d 

txM 2 M at A 



K o(2)<^ (7 



2 rff(»)w* 



a w 7 



iV-'iV 



M,7V — 1 N 



Dm,n, 



\ N l°M,iV ~ U U 

AT 



M,NJ f i 



where (f53"| has been used. The equation for 7 M follows by substituting (frJ2|) and (|55|) into (|54")) : 



.._E(1) 2 , 2 fl 2 S(l) 

7m ~ c 7m + c Pm7m 



C K (z) ^-V^ H ( 0)/ ^ _ j-ff(O)^*,* _ cjJgM\^'T/ H (°) 



»/-2 V'JV ' A" AT )^M,N Kf2 N ^M,N 

M - • /v i- 



N 



M 



N 



k (z) d w| 2 £(0 ) 2 B(0 v rB(o) rr * _ D $ X 

y-2A/"2 pit 2-^i I <Pn 1 N + t 11 IN J-^ALN J N l°M,jV U M,N) 



c 2 a 



a* A, 

^c 2 (9 



(pW-/>( )/to(«)rcoBfl)*jf#l 



M 2 M dt 



v 



(73) 



(74) 



Expressions for the other field components can be expressed in terms of "f M ^ and 7^- , their derivatives and lower 
order fields. From (JB5J), and ((Ml) one has 



V, 



,1/ 



E(l)l 

7m 



/i^W 2 



M JT> 



(p {1) -p (0) K (z)rcos6i)$M#l 



V, 



H(l) 



M 



a 



>P'/rv B(0) /? 2 _i_ ^^(o)'^ f-* _v B(0) r<7* — d* ^ 

p AT Pw + 7jv )^M,N V N l°M,iV U M,N) 

N 

tt (z) W g (Q) w ^ 1 

Tm + A/f 2 JV l°M,7V U M,N) h 



N 



T E(l) 



0: 



M 



M AT l 



M,N 



JV 

1 

^M, 



^ 1} <i>M#l 



T H(1) 



1 

4/ 



7m + 772- i (7w + "JV^N O-^iW 
■^M N L 



M.N 



T H{0) fr ,iy n * \ 
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3. General Solutions 

In WD II and W D 21 the problem of solving ([26]). ([27]) , (28]). |29|) has been reduced to solving an initial- value 

problem for the decoupled fields 7^°\ 7^^\ 7^°'' an d 7jv^- For some real constant cr > each satisfies a second- 
order hyperbolic partial differential equation in the independent variables (t, z), of the form: 

f-c 2 f" + c 2 a 2 f = g> (75) 
for some prescribed source function g. The causal solution of this partial differential equation for t > 0, with prescribed 



values of /(0, z) and /(0, z), has been exhaustively studied in the literature, see e.g. [12j. If the data and sources are 
sufficiently smooth the general solution may be expressed in the form 



where 



f(t,z) = H (T [f imt }(t,z)+l4g](t,z), (76) 



H a [f mlt ] (t,z):=Uf(0,z- ct) +f(0,z+ ct) 



2 



1 

2c 



z+ct 



d(mQJo(crV<?t 2 -(z-Q 2 ) 

z — ct 
z+ct 



* r^.w^-i-o 1 ). (77) 



and 



^^(t.z)^! [ Z+C(t ^ dtg(t',Q).h{cr^cm-t>y-(z-ty), 

Ac J0 Jz-c(t-t') 



(78) 



The functions f(0,z),f(0,z) constitute the initial t = Cauchy data in this solution and determine the H. a contri- 
bution above. Typically, in an accelerating device, lowest order contributions include externally applied piccewise 
established magnetostatic and RF fields that are together used to guide and accelerate charges along the beam tube. 
In the following we assume that all H a contributions to the field solutions arise in lowest order. 

E. Electromagnetic Power from Smooth Sources 

In the general situation it is seen that all zero and first order fields can be calculated in terms of finite range 
integrals involving Bcssel functions. It is of some interest to calculate how the instantaneous electromagnetic power 
flux depends on the first order curvature correction to that in a straight cylinder with smooth sources. This is obtained 
by integrating the Poynting vector field over the cross-section T> at an arbitrary point with coordinate z. In terms of 
the Poynting 2-form 

S (e, t, z, r, 9) := e (e, t, z, r, 9) A h (e, t, z, r, 6), (79) 

(2) (1) (1) 

such instantaneous power ^(e, £, z) is obtained by integrating S over T)\ 

(2) 



w(e,t,z):= / S (e,t,z,r,6) = / S (e,t,z,r,6), 
Jt> ( 2 > Jv ( 2 ) 

where S is the 2-form that does not contain dz in S . From the electromagnetic 1-forms, (1391) and (l40l) . it follows 

(2) (2) 

that 

w(e,t,z) 

= Y,^V$(e,t,z)I%,{e,t,z) f d$ N A#(dzAd$ N ,) 

N N' "' V 

+ ^^<( £ ,i,z)/f ( ,( £ ,M) / (#(dMd* M ))Arf$ M '. (80) 
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Taking into account the symmetries given by (|43|) and (|44| . and the fact that the sums over n',m' £ Z arc from — oo 
to oo, Eq. ([50]) can be re- written 



Writing 



one has 



and 



JV AT' 



x / d$ N A #(dz A (-I) - " -d*_„/ .,-(_„/) ) 

ZJzJ ^mC -1 )*" lH rn',p'(-rn') 



M AT 



x / (#(dzA(J$M))A(-l) m d*_ mV( _ m0 
= Y,Y. V nW> I d ®N A#(dzAd$^7) 

AT AT' *' :D 
A/ A/' 



w(e, i, z) = iuC) (t, z) + eu; (1 ) (t, z) + C(e 2 ), 
;(°)(t,z) = ^^^ (0) (t,z)4 (0) (t,z) / d^A^dS 

AT AT' J"D 

+ EE y M° ) (M)C ) (M) / (#<i*M)Ad%: 

A/ U> V 



» (1, (m) = ee(C(m)4 (0) (mkw 

AT JV' ^ 



x / r cos 9d§ N A #d$jy' 
+ (^ (0) (t,z)4 (1) (t,z) + V A f (1) (t,z)4 (0) (t,z)) / d$„A#d$ 



-t 'a/' 



+ EE{C ) (M)C ) (MW2) / rcos0(#d* M )Ad* 

M A/' ^ •^' D 

+ (V* (0) (*, ^ (*, *) + (*, ^ (*, *)) 

x J (#d* M ) Ad^M'j- 

With the aid of the symmetries, (|43p and (|33J) . and the fact that 

/ d $jv A #d $ap oc <J nV _ n , / (#d# M ) A« A/ ' oc <5 mV 



one obtains finally 



^ AT A/f J 



W Af 

where 3? represents the real part of its argument. 



(81) 



(82) 



(83) 
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Similarly the first order correction to the power is 



7« = ^|e^(^ (0) 4 (1) +^ (1) ^ (0) ) 

^ N 



/ J 0L M JVi M\ V M 1 M +V M 1 M 
M 



+«o( Z )EE4^; (o) 4 (o) 

JV N' 



M A/' ^ 

The analysis above is general and accommodates arbitrary smooth continuous conserved source currents. It should 
be stressed that in general one must specify how the prescribed sources should depend on both space and time variables 
so that their perturbative expansions in k can be determined. In the following sections we consider particular sources 
of relevance to the issues mentioned in the introduction. 

VI. SMOOTH LONGITUDINAL CONVECTIVE CURRENTS 

In this section a particular current source is considered. It is generated by an arbitrary smooth convective charge 
density p that is taken independent of k. It has components J r = Jg = 0, Jq(z — vt,r,9) = v p(z — vt,r,9) with v 

less than the speed of light. In this source model Jq ^ = Jq and all higher orders are taken zero. This is a contrived 
current source but serves as a comparison with the localised sources that will be used to model accelerated bunches 
of charge in subsequent sections. 

A. Lowest Order System with Smooth Longitudinal Currents 

The equations for 7^ and 7^ follow for these sources from above: 

-H(0) 2 H(Q)/I 2 2 -ff(O) n , Q r\ 



7^ (0) ~ c^T" + SPllT = - v 2 )p7- (86) 

where the conservation relation 

has been used and for k G Z>o we define 

:= / (87) 
Jv 

In terms of 7^°' ) and and the projected convective longitudinal sources the equations (|65|) . (|53|) . (|59|) and 

([55)) yield 

V E(0) _ 1 ( E(0)> 1 JQ) 

N ~ (3 N V N N N pY 2 

V H(0) _ M -H(Q) r „ Q v 

iT - -^(^ (0, +^). ( 9 °) 

^ (0) - ^7^ (0) ', (91) 
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and from (|83|) . (|88|) . (|89|) . (|90|) . and (|9T[) the lowest order contribution to the power flux becomes 



■w 



I (0) 12 --E(O) E(0)i 



N 

PN / m M 



B. First Order System with Smooth Longitudinal Currents 

The equations for 7 M and 7^ follow similarly. From |56|) . ([64]) . (j60|) one has 



~H(1) 2 , 2 2 

7m - c 7 M + c a M 7 M 



C 2 K (2)w 



1% _ V ) M ,N 772 W M 



A/f2 ^ V N )^M,N KA 2 X JV ^M,N 

JVl M N JVL M N 

/ H(0)// 2 T H(Q)/ ^ \ C 2 : k' q (z) ^-V / , H(0)r 
-ilN + a Nl N > E ALNj + ^ 

W N lT ) )Et I , N IT\CZ,N ~ D*^)}. (93) 

and dS2D, (ESD, HMD yield 

..£(1) 2 £(1)" , 2 fl 2 E(l) 

1m - c 1m + c PmIm 

~~ A/2 ^ W + M-'AT J^ALN \T2 2-^i N U M,N 

M N M N 



C 2 flKo{z 



N ^ 

C 2 K (z) ^ f 2 E{0} , g (Q)„, g, ^£(0)//^$ n $ 

775 2^ 1 yl J N V N + 7at )^M,N ~ V N l°M,iV ~ U M 

M N I 



<M JTJ ^ M 

+iT')Em,n - VT\C%, N ~ D^)}. (94) 

The right hand sides of these expressions for "f M ^ and 7 M can be rewritten in terms of the lower order fields 
1n^ ' In ^ anc ^ the sources. Thus the following terms in (|9"3"l) can be expressed as 

N Vn \p 2 n 1n kp n Pn r 

a N \c 

H(0), 2 rH{Q) _ „ H(0), 

In + a N 1 N — l 1n 
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Similarly the following terms in (194)) can be expressed as 



V H{0) ' + ui H{0) - 2-^ ff(0) ' 



'A' 



#2 rS(O) , y2»S(0) _ V (0)/ 
PAT- 4 TV + A" 7at - Jp P N ' 



AT 



fS(0) 
N 



AT 



A f 



#2 y B(0) 15(0)/ _ o -E(O)/ 

Pat + 7at — z 1n 



V, 



E(0) 



N 



N 



A/>Y 2 
To) 



(0) 

Pjv > 



Finally once 7^- and nave been determined from these decoupled equations it follows from (|66[) . (|56| . (|62|). 

(|64p that the remaining fields can be readily determined as: 



V, 



73(1) 



v, 



7/(1) 



JVf 



T E(1) 



jHQ) 



1 



I? 



73(1)/ 

7m 



M 

Ko{z)- 



pY 2 M 2 M Jv 



(p<°>/so(«)rcos0)$ w #l 



a/ 



A' 



/ .< | l JV Pw + 7/v )^M,N V N l°M,JV U M,N) 



a 



.HQ.) . Ko(z)x-*'irH(p) fri9 



M7M ■ M 2 



2~ JV l°A/.JV ^JM 



M.JVJ f' 



"pi 



M N I 



M.N 



T E(0) fr ,$ 



A 1 



JVf.JV ^M,N 



D 



A/ 



H(l)/ . k o(^)y^'// H(0)/ . 2 r - ff (°)\Z7* 

7m + 772-2^ 1 v/V + «jvV )^jvf,jv 
m m JV L 



T H(0) fr <S! n * \ 

J JV V°M,JV U M,NI 



Radiated Power Dependence on Local Curvature 



A straightforward but tedious calculation leads to the dependence on local curvature of the radiated power 



for 



the convective source model above, in terms of 7^°^ )7jv , 7/v j7jv > various overlap coefficients and projections 
of p. 

Each term in ([84|) can be reduced as follows: 



V- B 2 hi 2 V E{0) T E(1) - - V 

N N P N 



73(0)/ 

7jv 



N 2 yY 



(o) 
2 Pat 



,Y 2 ^ 



K (z) 



A/2 



Efl.N v (0) 



A/2 



°JV,M 



JV.M 



/3 2 



M 



V 2 • 73(0) . 

^ 7a/ + ITT Pa/ 



(0) 



A/ 



(95) 
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w-2 v E{l)jE(0) 



N 



N 

Af 2 



E(l)l 

In 



pY 2 Nf 



p (0) K (z)rcos6»$ A r#l 



AT Jv 



- E 



A I 
$ 

N.M 



2E 



N,AI 



r 1 ® _ 

U M,» M,N \ E(0)i 



P 2 



A I 



r<* n* \ (°) 

°iV,M U N,M \ Pm 



A I 



h! 2 



7m 



.£(0) 

In 



(o) 
v Pn 



(96) 



r H(Q) r H(l) 
' AI 1 M 



A I 



Ko( #E'(^ 



M 2 M H (0) J H{l)l 



M a M 



M 



A I 



N 



M.N 



7")* 

M,N U M,N \ H(0)' 



a 



N 



7jv 



(97) 



E« 2 

M 



A I 



Ml r V, 



AI V AI 



H(l) T H(Q) 



L M 



= »Y. M1, 



A I 



'AI 



■H(l) 

7m 



K o(z) st^' c m,n 



M AI 



E 



1J M,N 



N 



a 



■H(0) 

7jv 



JV 



H(0)l 

7m , 



(98) 



2__j N,N' v N 1 N' 

N N' 



Finally, 



AI AI' M M' m m ' 



VII. MOVING POINT CHARGE SOURCE 



In this section the above formalism is applied to the determination of the fields produced by the motion of a 
single point charge source moving in a curved beam pipe. A full account would involve generalising the formalism 
to accommodate fields from point distributional sources. Following established custom the formalism is extended 
here by modelling a point source as a moving Dirac distribution on U. One can then explicitly remove some of the 
integrations that arise in the smooth continuous charge source model. However a point charge implies singularities in 
the electromagnetic fields and these should not be ignored. The distributional charge density must model a charged 
particle with constant electric charge Qtot so we demand that 

p#i = Qtot- (ioi) 



u 

Furthermore the motion of the charge is maintained (e.g. by externally applied magnctostatic fields) on a curved path 
parallel to the design-orbit with curvature k(z) and constant speed v. These conditions are satisfied if 

p(e, t, z, x 1 ,x 2 ) = Q(e, tySixt - x 1 . )S(x 2 - x 2 . )S(z - vt), 

with 

Q(e, t) := %^ = Q tot + en Q {vt)x lfi Q tot + 0{e 2 ), (102) 

1 - 6K (vt)X lfi 
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in terms of the symbolic representation of the Cartesian three-dimensional Dirac distribution with moving point 
support at 

(0:1,0, X2,o,vt) = ( r o cos 80, ro sin#0j vt), determining the location of the point charge in IA at time t. 
Thus it follows from (fTU^]) that 

p = p(o) +e p« + 0(e 2 ), 

p (0) (z - Vt,X\,X2j = Qtot5{xi—Xifi)8(x2—X2fi)8(z-vt), 

p {1) (t 1 z 1 x 1 ,x 2 ) = Ko(vt)xiflQtotS(xi-x li0 )S(x 2 -X2,o)S(z-vt). 



In adapted coordinates 



(p#l)(e, t, z, r, 9) = Q(e, t)^ — ^5(6 - d Q )S{z - vt)rd r A d0. 

r 



The associated electric current components are J r = Jq = and 

Jo 0) = vp^ = vQ t otS{xi - x lfi )S(x2 - x 2 ,o)8{z - vt), (103) 

Jfi 1 ' = Vp (1) = VK (vt)xiflQ t ot 

xS(xx - xi fi )8(x 2 - x 2 ,o)5(z - vt). (104) 
Note that in this source model both Jq°\ Jq 1 '' and p(°\ p^ are non-zero with all higher orders zero. 

A. Lowest Order Fields from a Moving Point Charge 

Using this distributional model the projected sources in (|86|) can be evaluated as 

X ITz^ (2 " Vt > r ' V*"* 1 = QtotJn (^(«) 7) e- in9o 8'(z - vt). (105) 
so the solutions for 7^ and 7^ following from ([75]) are 

1 r\t,z)=H aN [iT )mit \(t^\ 
7 r ) (M)-^ JV [ 7 r iMt ](M) 

~^{c 2 ~v 2 )Q tot J n (*,(») ^) e- me "T N (t,z), 

where 

T N (t, z):= [ dt' f + i } dC6'(C - vt')J {(3 N ^c2(t-t'Y -(z-C) 2 ) 

and Ti aN is given by (|77|) . For i> > (See Fig[T]) it is shown in Appendix B that: 

• for z < —ct, (Outside the Regions Ri and R 2 ) 

Tjv(M)=0, (106) 

• for — ct < z < vt, (Region R 2 ) 

V (+ \ R f 1 '^ ^ (z-vtV^NS^zlt')) 

T N (t,z)=/3 N dt , 107) 

Jo s(t,z\t') 



for vt < z < ct, (Region Ri) 



T„(M) M» r'- ( ''^<'<^«M, (108) 
Jq s{t,z\t') 
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for z > ct, (Outside the regions Ri and R2) 



T N {t, z)=0, 



(109) 



where 



. . , t + z/c , , . t — z/c 

t' + (t,z) := — , , t'_(t,z)- 



1 + v/c 



1 — v/c' 



(110) 



and a(t, z\t') := ^c 2 {t - t') 2 - {z - vt'f 



The field components V^ (0) , v" {0) , I^ {0) , l" {0) follow from ®, flU, (71]), $M in terms of the solutions for 

ff(o) ^(0) 

IN ' 7w • 



"i 1 1 1 1 1 1 1 1 1 it — 1 r 



ct=-z 



ct=cz/v / / ct=z 





J I I I I I I I I I I I L 



FIG. 1: Domains for Tjv(t, z) and A a (t, z), Bcr(t, z) in flVII Bl Ri is the triangle in the upper right-hand corner (vt < z < ct). 
R2 is the adjacent large triangle {—ct < z < vt) where the source world-line has z — vt. 



1. Lowest Order Contribution to the Instantaneous Power from a Moving Point Charge 



We note from (|92[) that the lowest order contribution k/ ' to the instantaneous power contains terms derived from 
the projections 



p<f» {z - vt, r, 0)$at#1 = Q tot J n x q{n) -± e- mti °5(z - vt). 



'0 



(111) 



Such distributional contributions are absent in pipe cross-sections that do not contain the point source. Thus for 
sections with z ^ vt the instantaneous power is given to lowest order as 



N E(V)I.E{0) \ - Mjf .g(0) ff(0); 



M ~*M 

explicitly in terms of the above solutions. This is independent of the curvature re 



B. First Order Fields from a Moving Point Charge 



The explicit computation of the contribution of the radiated power to the next order, involving re, is somewhat 
more complicated since it requires the solutions to the zeroth-order fields as sources for the first-order fields as well 
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as evaluation of projected sources that now involve The equations for 7^^ and 7^ also contain additional 
source terms from these first order contributions to the sources. From ([55]) . (|64l) . ([ST))) one has 



2 h(i)" , 2 2 H(l) 



A//2 ^ N > 11 V N I M,N k A 2 2-j N M,N 

JVl M N JVl M N 



M M 



"(7JV + "nV ) E M,N j + I (7 



H(0)/ 



+ a jV^A^ ')-^M,iV — '(^M,iV — D M,iv)|) (112) 

and (32), (gS), dMJ) yield 

-E(l) 2 £(1)" , 2 fl 2 E(l) 

7m - c 7m + c /3m7m 

~ Kf2 2-~i ^ N ^ N > M,N Kf2 7. V JV "M,JV 

^M jv ^M jv 



fJ,C 2 K {z) W f 2 -^(0) y 2 »g(0)v F $ r-E(0)/ r $ n * 



A/" 2 

A « AT 



M.N , 



Kf 2 

JV M N 
C 2 



/ , 1 U% K JV +7jv ) M,N V N \ U M,N U M,N)( 

{p(D' -{pW K ' Q { z ) + p(V> KQ { z )) rC os6y¥^#l 



MC 2 ^7(TT f ,„ 2 ^£(0) e(0) V$ 

TFTPm H r?2 2^ 1 (Pn V N +7AT J^M.AT 

- ,V M • /V A/ N I 

The right hand sides of these expressions for 7 M and 7 M can be rewritten in terms of the lower order fields 
In ' In an( l the sources. Thus the following terms in (| 1 1 2|) can be expressed as 

T E{0)> _ llY 2 V E{Q) £*_V (0)/ + v 75)7 

J Af V- 1 V N ~ z l a2 7w + A/-2 o2 P/V 

\ HJV NP N 

,2 T VH(0) , r ff(0)/ _ 1 / 1 ,.H(0) , J?(0)// 



H(0)> , -,8 rff(0) _ o ff(0)/ 

7jv + ol n 1 n — z 7jv 
Similarly the following terms in (|113p can be expressed as 

V H(0), fH(0) _ 2 J*_ . H(0)/ 
"AT ^ M J AT — z 2 7iV 7 
"AT 



Pat-* at + A" 7w ~ TTTPn > 

n 



2 N — 02 \ t 11 IN K f2 "n 

PN \ • /v Ar 

PN v N ^ IN — Z 7at j^j-2 ^y2 < J N ' 



v E(oy _ 1 / E{Q)II _ (0)1 
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Finally once 7 M *' 1 '' an< ^ 7m nave been determined from these decoupled equations it follows from (|66l) . (|56[> . (|62p . 
([51)) that the remaining fields can be readily determined: 



V, 



E(l) 



M 



1 

w 



7m (1) ' 



V, 



H(l) 
M 



M 

«o(£) 
A/" 2 

■^M 

1 



(^-pW*so(«)rcosfl)* w #l 



U^iV "iV + 7 W J^M.iV - TV \ C M,N ~ U M,Nl 

AT ^ 



M 



— I M7m 



■ g(l) , W (^) W g(Q). ^ * I 

~/ / 'JV (,°M,JV — U M,N) p 



rB(l) 
'M 



Z M 



/3 2 



1 

M I 



/^ 2 7 M * 



, E (« (0) +^ (0) )^ 

7V M JV 



M,N 



T E(0) (r <s, 



V (1) 

AT2 PM 

■ /V M 



1 

T 

M L 



H(i), j_ Ko(z) s^'f f ..H(o)> ^ e^j^C ))^* 



7m 



■^M , 



E'{(7^ ( 



M,JV 



T H(0) (rr y _ n * v 

The distributional nature of the source permits evaluation of the integrals in (| 1 1 2[) and (|113p yielding 

- (V°v (z) +pW'«o(«)V«»«te#i 



= r cos 9 Q t otJ m (x p {m) —J e lm6 ° 

x S'(z — vt) I n {vt) — k (z) I — 5(z — vt)n' Q (z) 



and 



x f — K (wi)(5'(2 — ft) + KQ(ui)<5(z — ui) 



where «o(i»i) = dKo{z) / dz\ z — V f 
M 1 ) 



The solutions for j"/ 1 ' and 7^ follow similarly from (|76p . The source term in the equation for 7^ i '' is 



g$ 1 \t t z) = K$ 1) (t,z) + L$ 1 \t t z) t 



where, for kq varying with z, [l| 



#m (1) (M) 



M 



jV 



• -E(O)/ 
r M,N 02 7at 
Pat 



°M,AT ^M,N I 1 N H(Q)ll 



N 



7jv 



np* #(0)" 



Mir 



E 



m,at7jv 



iV 



2£ 



M,N 



°M,iV U M,N \ H(0)l 



7jv 



(114) 



Sit "* 2 ) 



2c 2 vQ tl 
M 2 M 



-K (Z)5 [Z - Vt) 2^ 772 m J " ( — ) e 



/ill),: 



c 2 vQ 



tot l 



M 2 M 



-K 



(vt)6(z - vt)J2 T7^r J « ( x 9Htt) 

N JVnPn v a ' 



-in&o 
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The source term in the equation for 7^ is 



where 





gf} 1] (t,z) = (t,z)+L E J 1 \t,z)+P M {1 \t,z) + (t,z) + SZ (1} (t,z), 



K 



M 



(t,z) 



Jy M 



C 2 K (z) ^-V f ~2fl . H (0)/ C M.N ~ D 



A/" 2 



N 



^-f^,w • 11 y 

u m,n")n 



<N 



c 2 [3 2 N 



M,N ..E(0) 

7JV 



2E 



M.N 



^M,N M,N \ E{Q)n 



N 



1 ^2 U M.n"?N 
a N 



/-<<E> 

°M,JV 



U M,N \ E(0)/ 



/? 2 



/zc 2 (c 2 + V 2 )Q t ot 



h! 2 



7jv 



Ko(z)<5'(z — ui) 



771$ 

J- 1 M.N 



°M,JV 



U M,N 



0N 



^m (1) (M) 



fic 2 v 2 Q tot 

aT 2 



r COS J m (%,(ro) 



x ^ K (wi)J(z — uf) — K (wt)(5'(z — ui) >, 



— ^2 Tq COS( 



A' 



M 



x - wt){K (wt) - k (z)} - 6(z - ut)K (z)^ , 



'm 



fJ-c 4 Q tt 
\T 2 

- ,V M 



M,N M,N 



M.N 







(115) 



In the next subsection. WII C[ the fields associated with an ultra-relativistic point source will be of interest. Then the 
contributions to the solution 7^°^ that depend on Qtot tend to zero and the solution 7^°^ depends only on external 
and static magnetic fields since there are no magnetic charges in existence. In this limit the source contributions 
K M ^ and K^ 1 ^ to g^ ^ and gf} drop out. 

In general, an analytic form for 7^ follows from ([77j) and (|78|) . by applying the integral operators H aM and T aM 

These source functions are simple functions of z, t multiplied 



to the source functions K^^ , £^ , Pjj , i^r , S^j 1 '. 
by (complex) numerical coefficients. Thus with 



h(t,z) 



Ko{z)S'(z — vt), 
K,' (vt)6(z — vt), 
Ko(vt)S'(z — vt), 



(116) 
(117) 
(118) 



we write the application of T a on them as 



Z*[/i](f,*) 



-A a (t,z)+B a {t,z), 

B a (t,z), 



where A a and B a , are calculated to be (See Fig[T|): 
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• for z < —ct, (Outside the regions R\ and R2) 

A*(t,z):=0, B a (t,z):=0, (119) 

• for — ct < z < vt, (Region R 2 ) 



I /•*+(*,*) 

A a {t, z):=— dt' K ' (vt')Jo(as(t, z\t')), 



2c 

(7 



B a (t,z) := - ^ dt ^—^ Mas(t,z\t )), 



(120) 



for vt < z < ct, (Region Ri) 

f t'_(M) 



Ar(t,z) -^Ycj dt'K' (vt')J (<7 S {t,z\t')), 

yj f \ CT f*'-^ , ,Ko(vt')(z-vt') T . , . ... 



(121) 



• for z > ct, (Outside the regions i?i and R2) 

A tT (t,z):=0, B a (t,z) :=0, (122) 

with t' + (t, z), t'_(t, z) and s(t,z\t') defined in ()110|) . Thus one sees that to 0(e), A a containes k' while B a containes 

K . 

The causal solution for 7^ generated by the source term above then follows from (|78|) and can be written: 

7 f f (1) (t,z) = l aM [9M (1) ](t,z) 

= l aM [KZ (1) ](t, z)+X au [L^](t, z)+X aM K {1) ](t, z). 
The first term on the right above depends on initial data but explicit expressions for the last two terms are: 

M M 
, G 



^mIV K^) = — V72~~ ( - A* M (M) + (M) 



N 

and 

„2 



^ot M \^M — KA 2 -^om (M) 



E'^^(^H?)e-^. (124) 

AT ^ NPN ^ " ' 



Similarly the causal solution for 7^- generated by the source term gf} can be written as: 

7m (1) (M) = l 0M [gZ {1) ](t,z) 

= l f3M [K*p]{t, z)+X Pm [L^](t, z)+X 0M [P^](t, z) 
+I 0M [R^](t,z)+I 0M [SZ {1) ], 

with 

„2/„2 1 „,2\ 



t ir E Wm \ A*c 2 (c 2 +v 2 )Q tot ( . , . , , 

X Pm l L M ](*, z ) = Jp - A Pm (*, z ) + %v/ (*> *0 



M 



E' 1 F $ °M,AT ~ M,N \ T ( r \ 
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T 0M [PM W ](t,z) = - Kr ^ tOt [Ap M (t,z)-B 0M (t,z) 



,-E(l) u+ ^ _ V£ 2 V 2 Qtot 
A/ 2 

xr cos 9 J m (xp{ m )— ) e~"™ e °, 
X /3 , / [i?f / (1) ](<,z) = 0, 

-r re^Wi^ > ^° A Qtot . ,, s 
X Pm\Sm KM) = —Tp. — ^ft/lM) 

E' 1 / p $ , Cf/.jV ~ D fl,N \ T ( r \ 



N 



E{ 1) 

The contribution Ip M [K M ]{t,z) depends on initial data and electric currents that vanish in the ultra-relativistic 
limit. 

1. First Order Contribution to the Instantaneous Power from a Moving Point Charge 

In terms of the lowest order modal solutions the instantaneous power iv^ 1 ' (t, z) for z ^ vt,t > can now be 
computed to 0(e 1 ) from 



V 2 \ " ^"n E{0)i(.E{1) , kq(z) ^> Cn,m D%,m -E(0)\ 
N Pn ^ N M Pm ' 



+/^ 2 e %{in (1)i + ^e' - c ^~ v )7rl7r 



, ., M M -H(0) \ H(l), «o(f) V^' ^p* ^MJV_^MJV\ i?(0)/\ 



EM 2 M ( . H (i) , k (z) v-^'Cm.w Df I N . m )\ h(o)/ 
-^2-( 7 M J + X^E^ -iN'jlAl 

M a M \ JVI M N a Af / 



E{0)/ . E(0) ■ H(0) H(Q)/ 

C. Ultra-relativistic Longitudinal Wake Potentials 

The wakefield formalism is designed to exploit the simplifications that arise by considering the unphysical (ultra- 
relativistic) limit obtained from charged sources moving at the speed of light. The resulting electromagnetic fields give 
rise to various wake-potcntials from which wake-impedances may be computed for ultra-relativistic charged bunches 
with prescribed charged distributions. The formalism is based on calculating the cmf induced on a spectator (test) 
ultra-relativistic point particle moving behind a leading ultra-relativistic charged particle with the same velocity but 
in general on a different orbit. Since the section above provides the electromagnetic fields for a point particle moving 
with arbitrary speed on an orbit (in general) off the tube axis (with transverse coordinates (ro,#o)) one may readily 
calculate the general longitudinal wake potential to the same order as the fields, by having the spectator charge, with 
transverse coordinates (r, 6), at a fixed longitudinal separation J > behind a right moving source particle. 

The definition [l3| of the ultra-relativistic longitudinal wake potential is taken as 

W^" e °\e,r,e,a):=-^- [°° dz £ L £±I, z , r , g) , (125) 



2T/2 
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where £z TO ' So He, t, z, r, 9) is the z-component of the electric field generated by the point source with speed v — c and 
charge Qtot ■ [201 ] Since the z-component of the total electric field is 



id z eU,t,z,r, 



N 



N 



one has 



with 



Thus, 



E^M (e, t, r, 0) = e £ 1^ (*, ^m(t, 9) + 0(e 2 ), 



M 



E(l 



From the orthogonality relation (|34|) . one calculates the projected longitudinal wake potentials: 



Hence 



(r- ,e ) 

Af 



W, 



II M 



(e,3) := / W,i ro ' eo) $M#l. 



A/" 2 



Eli) z + s 



Qtot J -SI 2 

From (|125|) . the definition of the ultra-relativistic longitudinal impedance is taken as 

dse l ^ /c w\{ ofio \e,r,9,I), 
and from the orthogonality relation (|34[) . one calculates the projected longitudinal impedances: 



Z^> 9 °\e,rAu) :=- 

r. 



(Z^' eo) ) M (e,u>):= I Z\^ aa >{e,r,e,w)^ M {r,6)m 



(126) 



(127) 



(128) 



To calculate (|128|) . one needs Tjy- (f, z). In the ultra-relativistic limit the expressions pi6p . dll7p . pi8j) become 



fi(t,z) 
h(t,z) 
h{t,z) 



Ko(z)6 (z — ct), 
K (ct)S(z — ct), 
Ko{ct)8 (z — ct). 



(129) 
(130) 
(131) 



One may write the application of the integral operator X a on these as 

l a [fi}(t,z) = -A a (t,z) + B a {t,z), 
lAh](t,z) = A a {t,z), 
l a [h](t,z) = B„(t,z). 

where A a and B a are given in the following domains: 
• for z < —ct, 

Aa(t,z):=0, B a (t, z) := 0, 



(132) 
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• for — ct < z < ct, 



A& (t , z 



t' c (t„ 



2c Jo 

G dt ^o(ct')(z - ct') 



• for z = ct, 



2c 



Ao(t,z) 



dt' K ' (ct')J (as c (t,z\t')), 
Ji{as c {t,z\t')), 



(133) 



4c 



z{t,z\t') 



dt'K' (ct')Jo(as c {t,z\t , )) 1 



B a (t,z) := f [* dt> K ° (ct 'l iZ ^ MasS,^')), 



(134) 



Ac J a s c {t,z\t>) 



for z > ct, 



A a (t,z) :=0, 



:=0, 



where 



and 



It immediately follows that 



s c (t,z\t') := ^c*{t-t'Y-{z-ct')\ 
t' c (t,z) :=i(i 
s c (t, ct|t') = 0. 



(135) 



(136) 



(137) 



The explicit form of A a (t, z) and B a (t, z) follows from a calculation similar to that outlined for A a {t, z) and B a (t, z) 
respectively in Appendix A. Thus A a and B a correspond to A a and B a respectively in the case when < v < c. 
Using the results in qVH Bl and (fT52]) . (fTM|) , (fTM|) , (fT^]) one finds 



7^ (1 V, z) = Z0 M + + RT + sZ {1) Kt, z), 



li r r 6o '(~A f3M +B PM )+p^r o, (A l3M -B, M ) + 



(r ,S ) 



'Af 



Af3 A 



j i(r<)fi a ) (rofio) (r ,8 )<. I , /.(ro,0o) (n> ,0 o) \ .5 



(138) 



where 



(r ,e ) 

M 



2^c 4 V-' 1 / p* 



A/2 



JV 



AI.N 



°M,iV ^ M.N 



Pm 



\xc 

A/2 



Qtot^o COS 0oJ m (x p( 



r 0\ -im9 



M — AT2 ^totZ^ A/2 ^ 



M,N 



°Af,iV 1J M,N 



P 2 



N 



(»«(n)^) 



-in9n 



(139) 
(140) 

(141) 



From these expressions one calculates (|128p and (|127[) . Thus 7^ and hence (|128[) can be expressed in terms of 
separate contributions from Kq(z) and k' {z). 
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1. Longitudinal Wake Potential for a Pipe with Piecewise Constant Curvature 

In the last section explicit formulae are given for the computation to leading order of the longitudinal wake potential 
in a pipe with arbitrary smooth curvature and |/c(,z)a| < 1. From such potentials one may calculate the longitudinal 
impedance to the same order. These expressions involve integrals of the curvature with Bcsscl functions and such 
integrals in general require numerical analysis. However in cases where segments of the beam pipe are connected by 
planar segments of arcs with constant radius of curvature (See FigJ5]) one may perform these integrals analytically 
and hence generate analytic expressions for the corresponding wake impedances. In principle there is an element of 
further approximation involved if one assumes that the tangent to the axial space-curve is discontinuous where the 
straight segment joins the curved segment. However bearing this in mind consider the case of an infinitely long planar 
pipe with axial curvature given by 

k (z) = (Q(z - z L ) - Q(z - z R )) k , 
where zl, zr, (0 < zl < zr), kq(^ 0) are constants and Q(z) is the Hcaviside function 

i 0, otherwise. (142) 
In this case one can calculate Aa,Bo- and write (|128|) in terms of known functions. 




FIG. 2: Profile of beam pipe with a segment of constant curvature. 

E( 1) 

One finds that the terms in the wake potential proportional to Ap M express the contributions to •y Af from the 
~ = Zl and z = zr and the terms proportional to Bp M express the contributions coming from the reg : "* 
vaturc is the constant kn. 



waiwiuiuiiLi a,u ^ — A, ^ o,±m A, — A, H ciiiii u±i^ u^iiiio piu^ui uiuiicu LJp M 

zl < z < zr where the curvature is the constant kq. 

The only non-zero contribution to A a {t, z) arises from the region 
integrand now follows from the relations 



ion 



transitions at z = zl and z = zr and the terms p 
zl < z < zr where the curvature is the constant kq. 

1 -ct < z < ct. The term involving K' (ct') in the 



and 



S(z/z ) = \z \5(z) yielding 



^- ^Q(z - z L ) - Q(z - ZR^j = 5(z - z L ) - S(z - z R ), 



A a (t, z) = J^ UZ) dt' [d (f -^)-8{t'- Z -f)^ J (as c (t, z\t')). 



Hence, see Fig EH 

• For (—ct < z < ct) n (zr < ct' c (t, z)), (Region R{) 

A*(t, z) = -^i J Q (<js c (t,z 



— j)-Jo((TS c (t,Z 



zr 
c 



for (-ct < z < ct) R (z R > ct' c (t, z)) n (z L < ct' c (t, z)), (Region R 2 ) 

Au(t, z) = ^ 



= S J [We [t,Z 



for (-ct < z < ct) n (zl > ct' c (t, z)), (Region R 3 ) 



A a (t,z) = 0. 
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FIG. 3: Integration ranges for A a (t, z) and B a (t, z). R\ is the triangle in the upper right-hand corner i.e.,(— ct < z < ct)C\(zn < 
ct' c (t,z)). R2 is the adjacent large rectangle: (— ct < z < ct) n (zr > ct' c (t,z)) n (zl < ct' c (t,z)), and R3 is the adjacent smaller 
rectangle: (—ct < z < ct) n (zl > c£' c (i, 2)). 



Since at z# = ct£, and zl = ct' c , 

s c I t, 2z R - ct 



zIL 

c 



0, and s c t, 2zl — ct 



zl 
c 



A simple calculation from l|134p . shows that with z = ct, 
• for zr < z, 



• for z R = z, 



for Z L < z < z R , 



• for z = Zl, 



• for z < Zl, 



Similarly since 



A a {t,z) = Q, 



A„(t,z) - 



■Act (t, z) ^ 3 i 



Acr(t, z) — ^ , 



Aa(t,z) = 0. 



B <7 (t,z) = o, 

for — ct < z < ct, and (z < zl) U (z > zr), we must evaluate 
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with — ct < z < ct and zl < z < zr where kq(z) ^= 0. 
Write the integral in B a as 



/•*'=(*'*) z-ct' 
K a {t,z):= I dt' z J^as^zlt')), 



and change variables, t' i— > s' c (t') with fixed i and 



s' c (t') := Vc 2 (t-i') 2 - (^-ct') s 



Then with 



Eq. (|143j) can be expressed 



d£ = ds 



c(z — Ct) ' 



2(z - ct) 



, (z - ct) 



2 „/2 



ds - 
- c 2c(^-ci) : 



-Ji(cr4), 



since 



s c (t,z\0) = V(ct) 2 ~z 2 , s c (t, z\t' c (t, z)) = 0. 
This integral can now be evaluated using 

/dzJi(az) = Jo(crz), I dzz 1 J\(az) = — J^iaz). 
a J a 

Hence 

2c{z-cty J^ [ct) 2_ z 2 

n 



K = r- / ds' c Ji(a4) 



l ("I) 



2c CT 

or with J o (0) = 1 and J 2 (0) = 0, 
A 



J (ct<) 



1 1 

2c(z - ct) 2 CT 



valid in the regime — ct < z < ct. 
Finally for z = ct, 



A„ = 0. 



Thus the function T^ 1 ' (t, z) is given in the indicated domains (See Fig[3] ) by: 

• for [z < —ct) U (z > ct), (Outside the regions R\, i?2 and R%) 

lM {1 \t,z) = Q, 

• for (-ct < z < ct) n (z L > ct' c (t, z)), (Region R 3 ) 

7^ (1) (M) = (4° A) -^ 0) )^A /3m (M), 

• for (-ct < z < ct) R (zx < ct^,(t, z)) n (zr > ct' c (t, z)), (Region i? 2 ) 
7^(M) = (-^)+^) +s ^))^J ( /3mSc ( M 

2c 



+(4 oA) -p£ 0i9o) )^?a /3m (m) j 



(143) 



(144) 



(145) 



(t, s) = ^- ( - 1 + J (CTV(ct) 2 -z 2 )) + 2 ffl ^ (tvW 1 ^), (146) 
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for (-ct < z < ct) n (zr < ct' c (t, z)), (Region i?i) 



x ^ Jo /3ms c t, z 



zr 
c 



+(4° A) -^ ' e ° ) )^A /3M (M). 



The explicit longitudinal wake potential in this case now follows from (|128|) and ()138l) : 



(147) 



Here, VVite^^ expresses the contributions from the abrupt transitions in curvature at z = zl and z — zr, while 
denotes that from the region of constant curvature zl < z < Zr. Splitting the range of integration in (|128D 



w 



(rofo) 
\\M,k 



according the domains associated with T^ 1 ^ (See FigfSJ) one has 

•A/? 



w (r ,6o) ( ~\ _ , % (j(ro,e ) _ (r ,6 ) _ Jr ,8 )\ 
" V \\M,edaes V ' ' ^ r~i \ M "m h M > 



|| M, edges 
zl-s/2 r z R +s/2 



s/2 



z L -s/2 



dz Af 



zr+s/2. 



,z)+0(e 2 ), 



(148) 



and 



W 



(r ,S ) 



M (Arofio) _ Jr Q ,6 )- 

r> \ m I'm 

Wtot 



dz Br-, 



,z I +0(e 2 ) 



(149) 




J I I L 



FIG. 4: Integration ranges for Wfy]l° d > gt , s . 



Again using the relation 



/ dzz^ 2 J v (V^) = ^=z (1+l/)/2 J„ +1 (V^), 

J V' 7 



(150) 
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with v — the integrals in (|148p are evaluated as 

pZL—s/2 />oo 
-s/2 



K 



ssl— s/2 Jzr+s/2 
R+s/2 

dzJ [ (3ms c 

Zl-s/2 



z + s 



z + s 



^ c JzR+1/2 



— ] ] - Jo ( /?a/s 



t) 



«0 



dz'J (fi M V^z f ) 



-Jx(V2(3 M Z), 



V2c 2 Pm 

independent of Z£ and zr. The integral involving B a in (I149[) can be similarly evaluated using (|150[) with ^ = 2 



(151) 



2c 



4c 2 



2c 



dzKp h 



z + s 



2c/3 
2ZS + S 2 

' 2C&MS 1 

Pmvs 



M 



-1 + J (PmV^z+¥) 



J 2 (/3a/V2sz + S 2 ) 



O z '3/2 , 

+ ^^J 3 (/?MV2Sy) 



-| z'=zr+s/2 
z'=zl+s/2 



Finally from ([T3g]) . (fT5"Tj) one has 



(ro,eo) 



K n(ro,9 ) ,*(n>,0o) "(r ,e ) 



M,edges V 6 ' s ) ~ e 



)Ji(y2/3 A/ 5) + 0(e 2 ), 



and from (fHgj) , (fT52|l 



KV ||M,k [ - t ' b > ~ 6 4 V'M % J 



(zr - Zi) + 



72 
PmVs 



x ^ \/ z fl + - Ji /3m 4 / 2s Zfi + - 



+ 



2 (** + f ) 



Js Pm\ 2s z R + 



ZL + 77 Ji ft/W2s zl + 



2 zl 



,3/2 



-J 3 /?mW2s z l + 



+ 0(e 2 ), 



(152) 
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where, from (I139|) . (I140|) . (|141[) wc introduce the abbreviations 



'(r ,e ) 

M 



Nj[ j(r ,e ) 

c*Q to t M 



N 



N 



p $ ~MjN M.N 

Pn 



T ( r °\ 

x j n yx q ( n ) —J e 



r \ m6o 



Pm 



b M 



^M Jr ,B ) 



C 2 Qtot 



•Nil >o,flo) 

c 2 Qtot 1 



1 



= ^Y, TfT ( - E M, 



r ^ -imB a 

a 



AT 



°M,iV U M.N 



T ( r \ 



f 0\ -in8 



It is worth noting that the expressions for the wake potentials are independent of Qtot- 
With the following dimensionless variables for some length L 

Ko := Lko, s :— — , p M := LpM, Zi? := ~F/' Zl '' = ~L' 
one may introduce the dimensionless quantities 



/3m 



Cm,2(s) := ko 



- zl) 



V2 

PmVs^ 



x i \/2k+ o J i Pm\ Zs\ zr + 



2 + f ) 



3/2 



^3 @M\ Zs[ Zr 



Zl + -Ji[ Pm\ /2s zl + - 



2(^ + f) 



3/2 



-Ja I Pm\1'Zs\ Z L + - 



in terms of which 



>V 



(r ,e ) 

||M,edges 



( e >s) = — Cm,i(s)(^ m 



(r o ,0o) _ ^('"o.So) _ g(^o,6o)\ 



Pm 



W 



(r ,8o) 
\M,k 



(e,s) = -XmA^)(1 ( m° M -Pm qA)) ) 



Natural choices for L include L = a or L = zr — zl. In Fig|5]we plot Cm,i and Cm, 2 for the choice 

k = 1, /?a/ = 1, zr = 2, z L = 1. 
In the regime s"^> 1, Cm, 2 tends to kq(zr — %£)■ 



VIII. SMOOTH CONVECTED LOCALIZED BUNCHES WITH FIXED TOTAL CHARGE 



In the previous section attention was concentrated on the ultrarelativistic limit. In this section, by contrast with 
WI[ wc construct a moving source model with finite total charge Qtot: and a smooth charge density satisfying (I101|) 
moving at less than the speed of light. Large numbers of charged particles moving with a common axial velocity may 
be modelled by a localised smooth distribution of electric charge with a prescribed convective axial velocity field with 
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constant longitudinal speed v < Cq (independent of the local curvature k(z) ), charge density p{z — vt, r, 8) and current 
density components J r = Jg = 0, Jq{z — vt,r,8) = v p(z — vt,r,9). The localized charge density profile is thereby 
maintained as a function of the arc-length parameter z as it travels along the beam pipe. Note that in this source 
model Jq ^ = Jo and all higher orders are taken zero. 
We assume here that p can be expressed as 

p{e,t, z,X!,x 2 ) = Q(e,t)p ± (xi,x 2 )p"(z - vt), 

where v is given (v < c), p ± (x\, x 2 ) and p^(z — vt) are arbitrary smooth functions subject to 

/oo 
dzpHz-vt) = 1, 
-oo 

and 

Q(e,t) ~ .. Qto < p— = Qtot + eQ tot { P K a ) z (t){x lP ^) v + 0(e 2 ), 

1 - e{p"K ) z (t)\ x iP )v 



with 



/oo 
dz pW(z - vt)no(z), 
OO 

(x 1 p J -)- D := / x 1 p- L (x 1 ,x 2 )#l- 



Correspondingly, we define 

p(°\z - vt,xi,x 2 ) := QtotP ± (xi,x 2 )p^(z - vt), 

p (1) (t, z, X!,x 2 ) := Qtot (/»" K ) z (t) (xt p x ) v p X x 2 )p" (z - vt). 

The currents and are defined by 

j{f\z — vt,r,8) := vp^°\z — vt, x\,x 2 ), 
J$\t, z,r, 6) := vp {1 \t, z,x x ,x 2 ). 
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A. Lowest Order Fields from Moving Smooth Convected Localized Charge 

The equations for 7^ and 7jy are given by (|85|) and (|86|) . The causal solutions are given by 



7JV (M) = ?W7w ](^2)-T72-( C ~ W ) J PAPN KM), 

• /V 7V 

where is given by (|87|) and 

p { ° ) '(t,z) = Q tot ptt'(z-vt) [ p ± (x 1 ,x 2 )^ N #l. (153) 



In the ultra-relativistic limit, v — > c, the second term in (t, z) tends to zero. 

B. First Order Fields from Moving Convected Localized Charge 

The computation of the longitudinal fields follows along the lines detailed above for the point source. The distri- 
butional source is simply replaced by the smooth sources (|153|) at each respective order. For a general such source 
one obtains a system of fully coupled modal equations. 

The causal solution to j M for 7 M can be written as 

7 f 7 (1) (t, z) = l aM [Kl {X) ]{t, z) + l au (t, z). 

^H(l) U x „. . . mc nf H(0) , E(0) ■ o _ oo , ^H(l) 



As before the functional form of K M { '(t, z) given in terms of y N and y N is the same as for K AI '{t, z) defined 

>H(Y 
J M 



in (|114p but now, L M (t,z) is calculated to be 



L M (t, z) := 2- rT Ka{z)2_^ . , 2 r2 p n {t, z). 



Similarly the causal solution for 7 M is written 



7 f/ 1) (i, z) = 1 i3m [KfV>]{t t z)+l l3M (L^](t, z)+1 3m [P* (1 >](t, z) 

+1 Pm [R E M {1) ](t, z) + 1 3m [S^] (t, z) . (154) 

with the functional form of (t, z) written in terms of 7 jv^ and 7^°^ analogous to that of (t, z) given in 

(fTT5)) and 



>E(1),. v pC 2 (c 2 +V 2 ) W 1 

J M (M) ■= 772 K o(2)2^T? 

JV M N JV N 

M AT >-> 



2 



. ^M,N ^M.N (0)1,, \ 

x i - Em,n + —p — }Pn (*> z )> 



• /V A/ 



where 



■ /V M N ^ P N 



Pm = Qtot(xiP^)v^- t ^(p ll K Q ) z (t)pU(z - wt) j jT ^(xx.aJaJ^Af^l, 
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and 



p 



tot{p^Ko) z (t){xip ± )vp^ \z - vt) j p ± (xi,x 2 )®M#l 

Jv 

Qtot\no{z)p ll {z -vt) + K {z)p ll '{z -vt) \ / x 1 p ± (x 1) x 2 )$M#1- 



1. Axially Symmetric Smooth Convected Localized Charge Distribution 

If the transverse distribution depends only on r, expressions for the electromagnetic 1-forms simplify. The source 
is axially symmetric if 

p ± (x u x 2 )=K(r), (155) 

where IZ(r) is a smooth function satisfying 

/ drrTl(r) = —. (156) 
Jo 27r 

It follows that 

pP'(t, z) = 2TTQ tot S nfi p^(z - vt) J drrTl(r)J (^(o)") • (157) 

The expressions for 7^ and 7^ 1 ' simplify and explicit expressions for the overlap coefficients can be computed in 
terms of transverse projections of the radial profile lZ(r). We consider this case in simplifying the situation. From 
(|157p . it follows immediately that 

ifr (1) (t, z) = +<S m ,-i) ™L V Q to tKo(z)p^'(z - vt) 



\f 2 

- ,V M 



x K f2 M '°«2 (0) / drrn{r)J U q[0) r -), 

• /v n nimPn „(n\ JO v a/ 



<j(0)eN 



and 



?E(X) U \ fx , X ,^PC 2 (c 2 +V 2 )Q tot i \ II// ,x 

L M ' '(t,z) = (d ro ,i +d m ,-i) Ko(z)/>" (z-vi) 

^ M 

x x ■ , _ E <S> , °M,0,q(0) ^Af,0 )3 (0) 



XI ^2 

x / drrK(r)J (z<?(0)^) 



Kf2 } ^M,0,?(0) 1" " 02 
9 (0)SN 0,g(0) ^ ^0, g (0) 



Then, since 



one has 



and therefore, 



(V)b = 0, (158) 



p$(M) = 0. 



Pf/ 1) (i,z) = 0. 
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Finally, using (| 1 58[) and 

Xi/9- l $m#1 = T(5 m i + 5 m i) / drr 2 TZ(r)J m (x p i m \- J , 



one has 



- (V°V (z) + P (0) '«o(«))asi}^7#l 
x y drr 2 TZ(r)J m (^(m)^) • 



It follows that 



i? M ^ '{t, z) = -j-f2-QtotT(6m,-l + <W) 



A/" 2 



x ^K (^)p"'(« - Vt) + K (z)p"(z - vt)j J drr 2 K(r)J m (®p(m)~^ 



and 



V- / $ °M,0,g(0) ^M,0,<?(0) 

X 2^ 772 \ ~ fi «A?(0) H fl2 
g(0)£N °.9(°) V M °.9(°) 

In summary, for a smooth axially symmetric charge distribution, 7 M V and 7 M can be written as 

% z) = l aM [A* (1) ](t, z) + (S m ,-t + S m<1 )l 0M [i£ (1) ](t, z), 

and 

+(*m,-i + <W)^m [2m (1) + R E m ] + S^Kt, z), 



using 

(<5m,-l + <5m : l) 2 = (<5rn,-l + #m,l)- 

IX. CONCLUSIONS 

This paper offers an analytic perturbative approach to the computation of electromagnetic fields generated by a 
variety of charged sources moving with prescribed motions in a perfectly conducting beam pipe of radius a with 
planar curvature k(z). Results are given in terms of expressions involving powers of |a«;(z)| <g; 1 and |a 2 K'(z)|. It 
has included a discussion of ultra-relativistic longitudinal wake potentials from which pipe impedances induced by 
k(z) 7^ can be calculated. This has been explicitly illustrated for pipes with piecewise constant curvature modelling 
pipes with straight segments linked by circular arcs of (arbitrary) finite length. 

There are a number of extensions that follow from this work. They include the effects of non-planarity and non- 
circular cross-sections, both varying with length along the pipe, resistive pipe boundary conditions, generalisations to 
dielectric channels and the computation of fields satisfying periodic boundary conditions in cyclic machines. Although 
more challenging the use of the geometric perturbation technique presented here is immediately applicable to these 
problems and will be discussed elsewhere. 
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Appendix A 



From the expressions for the orthonormal coframes (|30[) in adapted coordinates it is straightforward to derive the 
decompositions: 

d#(dz Ad# M ) = -a^(l + ert O (z)rcos0)*M#l 
-£K Q (z)(#dty M ) A (cos 9dr - rsinddO) 
—eK' (z)r cos 

d#(dz A d$Ar) = -/3^(1 + eK (z)rCOs6>)<I>jV#l 
— eKo(z)(#<i$7v) A (cosOdr - rsin9d9) 
-eK' (z)r cos 6 (#d$ N ) Adz + C(e 2 ), 

d#d^ M = -(I - £Ko(z)r cos 9)a 2 M ^ M (#l) A dz 

-eK (z)(^d^M) A (cosOdz l\dr — r sin 9d z A d 9) , 

d#d$Ar = -(1 - £K (2 : )f , cos6l)/? 2 ^$Ar(#l) A dz 

— eK (-z)(#d ( E > 7v) A (cos9dz A dr - rsin9dz A d6). 

Furthermore if ^(r, 0) represents either $Ar(r, 0) or ^ m(t, 9) one has 

#d ip = (dMrd9- -(de^dr, 
r 

#dip = (1 - EK (z)r cos 9) (#dip) Adz, 



#(dz Adip) = 



1 — eKo(,z)r cos ' 



and 



In terms of the symbol 



— ife 



#dz = - #1. 

1 — eK (z)rcos# 



(J 7 , Q,x,y):= J drr k T m [x p{m) ~) 0„ ^ (n) ~ 



the following overlap coefficients are defined: 

2-7T 



7r((5„_ m+ i,o + ^n-m-i.o)-^- / drrJ m (zp (m) -) ^ (^(n)-) 



o 

a;', 

= 1"((Jn,m— 1 + <5ri,?ri+l) ~ ^M,N\^'^ ,X,x). 



and 
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Df IN := [ *Af(^*iv)#l = 7rn(5„ >m _i - 5 n , m+ i)E° MN (J, J,x',x'), 
Jv r 

e mn '■= / rcos0*Af*iv#l = 7r(5„, m _i + 5 n . rn+ i)S 2 M N (J, J, x',x'), 
Jv 

F MN '■= I rCOSOd^M A #d*JV = 7r(<5„, m _i + £n,m+l) 

x[ - - — ^, MJ y(J ,J , x , x ) + mnzP M N ( J, J, x , x ) ) , 
\ a a / 

Cmn '■= / cos(9¥^(9 r $iv)#l = 7r(5„, m _i + <5„, m+ i)-^-SL iV (J, J',x,x), 



v a 

-°m w : = [ ^^$M(de$N)#l = Trn(d n . m -i - 5 n . m+ i)E° M N (J, J,x,x), 
Jv r 

e mn '■= / rcos0$ M $jv#l = 7r(5„ >m _i + 5„.„ l+ i)S^ ^(J, J,x,x), 
Jv 

F M N '■= / rCOsQd^M A#d® N = ^(S n ,m,-l +S n ,m+l) 

Jv 

x z, M N (J , J , x, x) + mn^ M . N (J, J, x, x) 

\ a a 



G?;*„ := / rcosOd^M A <i$ 



Wr(<5n,m-1 + <$n,m+l) ( S M,Jv(-^'j ^ x 'i s ) 



-TO 9 ^ - aLjf(J, </', a;', x) I , 
a 



:= / rcos6d$M A 



Wr(<5 n , m _l + £n,m+l)^1 V{ ^ l) ^} MN {J' , J, X, x') 



+TO 3™ pj(J, J ,X,x) 

a 



Appendix B 

In this appendix the derivation of the expressions for T^(t, z) in (|106j) . ()10T|) . (|108|) . (|109j) . and A a (t, z), B a it-, z) 
in 1(115]) . IPD]) . (I2T1) , (PH is outlined. Since 

dzf(z)6'{z - zq) = -f'(zo), (159) 



dz/(z)<5(z-z ) = /(z ), (160) 



for a smooth function / and zq contained within Eel, double integrals involving smooth functions and 6 or 5' can 
be reduced to a single integral. For T^r, A a and B a , the triangular integration range defined in the (t',() plane by 
< t' < t and z — c(t — t') < ( < z + c(t — t') is reduced to a line segment specified by the support, t' = £/«, of the 
distribution S(( — vt). The slope of this line can only vary between ^ and infinity. The value of this slope determines 
which side of the bounding triangle the line intersects: see Fig[6ja) and Fig[6jb). For v > if z > vt the line segment 
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intersects the line t' = (t — -) + - at (£ = -f-^, t' = ^z^)- If z < vt the line segment intersects the line t' = (t + 
at (C = %P£,t = 2±M). When either 

z < —ct, or 2 > ct, 

there is no intersection of the line t' = (/v with the triangle and so Tw = yielding (I106|) and (|109p . When 

— ct < z < ct, 

there is an intersection of the line segment and the triangle, as shown in Fig[6]Ja), since 1/v > 1/c. 
If 

K(t, z)~ f dt' f + i } dCS'(C ~ vt)f(t, z\t', C), 



c ' c 



and 



it follows that if — ct < z < vt 



and 



where 



Similarly if vt < z < ct 



and 



where 



J(t,z) := / dt! 



:+c(i-t') 



d{6((-vt)f(t,z\t',0, 



J z-c(t-t') 



J 



t'(t,z) 



dt f(t, z\t', vt' 



K 



{t,z\t',0, 



C=vt' 



I , N t + Z/C 

t' + (t,z) := — -L-. 

1 + v/c 



J = 



t'_(t,z) 



dt'f(t,z\t',vt') 



K 



L 



dt dC 



f{t,z\t',(), 



C=vt' 



t'_(t,z) 



t — z/c 
1 — v/c 



(161) 



(162) 



Using the Bessel relation Jq(z) = —J\(z) the integrals p07p and p08[) follow immediately from p61[) and (|162[) . 
Similarly the integral representations of A a (t, z) and B a {t, z) follow from those for K. and J . 
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